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Abstract 



In this article, we discuss a couple of nonlinear Galerkin method (NLG) in finite clement 
set up for viscoelastic fluid flow, mainly equations of motion arising in the flow of 2D Oldroyd 
model. We obtain improved error estimate in L°°(L 2 ) norm, which is optimal in nature, for 
linear finite element approximation, in view of the error estimate available in literature, in 
L 2 (H 1 ) norm. 
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1 Introduction 

Nonlinear Galerkin (NLG) method has been introduced by Marion and Temam in [15J as a 
means of turbulence modeling, in order to handle long-time integration of dissipative evolution 
partial differential equations. In a classical Galerkin method (CGM), the equation is projected 
to some finite dimensional space, thereby neglecting the orthogonal space. But it is well-known 
that nonlinear dynamics is sensitive to initial data; a small change in the initial data, after a 
large time, may result in a significant change in the system. Hence, for large time integration, 
it is only appropriate to capture the effect of some the neglected terms. The idea behind NLG 
is to take into account the effect of the neglected terms in the long run. 

The method in |15j is based on the eigenvectors of the underlying linear elliptic operator. 
Later in |16| , it is expanded to more general bases and more specifically to finite elements. Unlike 
Galerkin finite element method, where we work with a finite element space Jh, with grid size 
h, in nonlinear Galerkin finite element method, we work with two finite element spaces, namely 
Jh on a coarser grid and jjf , a suitable complement of Jh in Jh, respectively. It corresponds 
to spatial splitting of the unknown u € Jh as 



In other words, we work with two set of equations in nonlinear Galerkin, one in y and the other 
in z. But the equation in z, which is solved in a finer grid, is generally simplified (linearized and 
time derivative is discarded). And the full equation in y is solved only in a coarser grid. In a 
way, nonlinear Galerkin method, in the context of finite element, results in a two-level scheme. 

In this article, we consider semi-discrete nonlinear Galerkin approximations to the following 
system of equations of motion arising in the Oldroyd model (see (|5j)) of order one: 



u = y + zeJ H + J^. 




du 



Vu-fiAu- / fi(t - t)Au(x,t) dr + Vp = i{x,t) 



x G O, t > 




o 



with incompressibility condition 



(1.2) 



V ■ u = 



ie!l, t > o 



and initial and boundary conditions 

(1.3) u(x,0)=uo in ft, u = 0, on dfl, t > 0. 

Here, Q is a bounded domain in R 2 with boundary dfl, fi = 2k\~ 1 > and the kernel (3(t) = 
7exp(— St), where 7 = 2\~ l (v — kX^ 1 ) > and S = A -1 > 0. For further details, we would like 
to refer to [5] and references therein. 

We are interested in optimal error estimate of the nonlinear Galerkin approximations of the 
velocity of (jl.lj) - (jl.3p . Note that the results obtained here for the above problem are also valid 
for Navier-Stokes' equations. Error estimations of nonlinear Galerkin methods in mixed finite 
element set up for Navier-Stokes' equations are carried out in pQ. Both nonlinearity and time 
dependence are treated on coarse space. It is proved that 

\\(u h -u h )(t)\\ m{Q) <c(t)H 2 , 
\\(p h -Ph)(t)h H n)<c(t)H 2 , 

where (\ih,Ph) is the classical Galerkin approximation and (u h ,p h ) is the nonlinear Galerkin 
approximation. These results are improved to 0(i? 3 ) in |17| . and also similar result is obtained 
for L°°(L 2 )-norm error estimate for velocity approximation, but only for semi-linear parabolic 
problems. The author feels that it is not straight forward to carry forward these results to 
Navier-Stokes equations or to the equations being considered here. The proof of [17] depends 
on one important estimate involving the nonlinear term /(u), namely 

K/(y/i + z ft) - f(yh),x)\ < H z fellL2(n)llxllL2(c)- 

(For the notation, kindly refer to [17] ) But it may not be possible to establish similar estimate 
for the nonlinear term /(u) = (u • V)u. On the other hand, it is observed that the improvement 
in the order of convergence is due to the fact that the splitting in space is done on the basis of 
L 2 projections, unlike previous approaches, where the splitting is based on hierarchical basis. 
Similar approach is adopted in this article. 

No further improvement is observed by the author for (piecewise) linear finite element dis- 
cretization and with forcing term f € L 2 ($7). Similar results are observed in p 7 , 8j, but for 
f € H 1 (J7) and for Navier-Stokes equations. In [18J, various nonlinear Galerkin finite elements 
are studied in depth for one-dimensional problems and similar results are obtained, although 
termed as optimal in nature. For example, the method is studied for piecewise polynomials of 
degree 2n — 1. And the error estimates obtained in energy-norm and L 2 -norm are as follows 

\\(u h - u h )(t)\\ L2{Q) < c^mmir^H^- 1 ,^^ 2 - 1 }, 

respectively. For n = 1, we have the piecewise linear finite element approximation and error 
estimates are of 0{H 2 ) and 0(H 5 / 2 ), in energy-norm and L -norm, respectively. Later on 
He et. al have studied NLG and modified NLG in both spectral and finite element set ups, 
[3j [9j QUI El [12] to name a few. In [3], (modified) spectral nonlinear Galerkin method is applied 
for the problem fjl . 1 [) - (|1 .3j) with periodic boundary condition. And in [101 1131 [12] , convergence 
and stability is analyzed for the fully discrete NLG. 

Recently in p3], a new projection is employed for a two-level finite element for Navier- 
Stokes equations. Error estimates of 0((log h)^' 2 H 3 ) and 0((\ogh) l l 2 H 4 ), in energy-norm 
and L 2 -norm, respectively, are obtained. The finite dimensional case of the Brezis-Gallouet 
inequality plays a crucial role apart from the projection in the error analysis and also the rea- 
son for the appearance of the logarithmic term. Note that here the forcing term is taken in 
L°°(0,T; (L 2 (fi)) 2 ) n L 2 (0,T;Hj(f})). The projection in [H] is based on fully discretization; 
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finite element for space discretization and backward Euler for time discretization. However, 
Crank-Nicolson scheme is used for the problem. 

It is the observation of the author that although the energy-norm error estimate for NLG is 
of order 0(-ff 3 ), there is no L 2 -norm error estimate of order 0(H 4 ). As is mentioned in [T7], it 
is an open question and the author tries to resolve this for a couple of (modified) NLG methods 
in this article. Although, this improved error estimate is proved for the problem (|l.ip - (|1.3p . it 
will go through for Navier-Stokes equations. 

The author would like to remark here that two NLG methods are considered here only for 
theoretical purpose. In the first NLG, nonlinearity is preserved in both the equations and as 
a result, is not physically viable due to very high cost of implementation. But by comparing 
the error analyses for both the methods, we can identify the key factors for the gain in order 
of convergence. In both the cases nonlinearity (in error) is separated and treated only after the 
linearized error is estimated. It is observed that the nonlinear terms play an important role, as 
they can limit the order of convergence, see Remark 15.21 

It is to be noted that the approach here is an heuristic one and the author does not claim 
that NLG performs better than classical Galerkin. The arguments in this regard are settled to 
certain extent in the articles \TE[ |6], to name a few. 

The article is organized as follows. In section 2, we briefly review the Oldroyd model. Section 
3 deals with the classical Galerkin approximation. In Section 4, we present the nonlinear Galerkin 
modelsand finally in Section 5, we discuss the error analysis. 

2 Preliminaries 

For our subsequent use, we denote by bold face letters the M 2 -valued function space such as 

H* = [2#(n)] 3 , L 2 = [LW. 
Note that Hj is equipped with a norm 

2 2 

l|Vv|| = (J2 (djV t ,d jVl )) 1/2 = {^2(Vv u V Vl )) 1/2 . 

i,j=l i=l 

Further, we introduce divergence free function spaces: 

Ji = {4> £ Hq : V • = 0} 
J = {4> € L 2 : V • <f> = in f2, • ti\qq = holds weakly}, 

where n is the outward normal to the boundary d£l and <fi ■ u\qq = should be understood in 
the sense of trace in H" 1 / 2 ^), see [19]. For any Banach space X, let L p (0,T; X) denote the 
space of measurable X -valued functions cj> on (0, T) such that 

f \\(f>(t)\\ p x dt < oo ifl<p<oo, 
Jo 

and for p = oo 

ess sup \\4>(t)\\x < oo if p = oo. 

0<t<T 

Through out this paper, we make the following assumptions: 

(Al). For g € L 2 , let the unique pair of solutions {v G J i,q £ L 2 /R} for the steady state 
Stokes problem 

-Av + Vg = g, V-v = inn, v| an = 0, 

3 



satisfy the following regularity result 

!I V I|2 + h\\m/R < C||g||- 

(A2). The initial velocity uq and the external force f satisfy for positive constant Mq, and for 
T with < T < oo 

u G Ji, f,f t a°°(0,T;L 2 ) with ||u ||i < M , sup {||f||, ||f t ||} < M . 

0<t<T 

Before going into details, let us introduce weak formulation of (|l.ip - (|1.3p . Find u(i) G Ji, t > 
such that 

(2.1) (u t , 4>) + /i(Vu, V0) + (u • Vu, 4>) + f p{t - s)(Vu(s), Vd>)ds = (f, 0), V0 G Ji. 

JO 

We would like to note the positivity property of the integral. For a proof, we would like to refer 

to [am. 

Lemma 2.1. For arbitrary a > 0, i* > and </> G L 2 (0,i*) ; £/ie following positive definite 
property holds 

J (^j exp[-a(i- s)](f>(s)ds\ </>(t)dt>0. 

We will use this non-negativity property frequently in this article without exclusive mention 
of the Lemma. Further, for existence and uniqueness and the regularity of the solution of the 
problem (|2.ip . we refer to [51 0] and references cited therein. 

3 Classical Galerkin Method 

From now on, we denote h with < h < 1 to be a real positive discretization parameter tending 
to zero. Let and L^, < h < 1 be two family of finite dimensional subspaces of Hq and 
L 2 , respectively, approximating velocity vector and the pressure. Assume that the following 
approximation properties are satisfied for the spaces and L^: 

(Bl) For each w G Hq n H 2 and q G H 1 /R there exist approximations ihW G and jhQ G 
such that 

||w - i h w\\ + h\\V(w - ih,w)\\ < i^o/i 2 ||w|| 2 , \\q - jh,q\\ < -^o^lklli- 
Further, suppose that the following inverse hypothesis holds for w/j G H^: 
(3.1) HVwfcll < Koh^Wwhl 

For defining the Galerkin approximations, set for v,w,0 G Hq, 

o(v,0) = (Vv,V0) 

and 

1 1 

b(y, w, <t>) = -(v • Vw, d>) - -(v • V0, w). 

Note that the operator b(-, ■, ■) preserves the antisymmetric property of the original nonlinear 
term, that is, 

b(yh, ™h, w h ) = \/v h ,-w h eH h . 

In order to consider a discrete space, analogous to Ji, we impose the discrete incompressibility 
condition on H^ and call it as J^. Thus, we define Jh, as 

3h = {vh G H h ■ (Xh, V • v h ) = Mxh G L h }. 
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Note that Jh is not a subspace of Ji. With Jh as above, we now introduce the following weak 
formulation as: Find Uft(i) G Jh such that u/^O) = uo^ and for cf) h G Jh,t > 

(3.2) (u ht ,4> h ) + fia(u h ,4> h ) + [ Pit - s)a(u h (s),4> h ) ds = -b(u h ,u h ,(f> h ) + (f,<f> h ). 

J o 

Since Jh is finite dimensional, the problem (13. 2D leads to a system of nonlinear integro-differential 
equations. For global existence of a unique solution of (|3.2p . we refer to OH]. 
Moreover, we also assume that the following approximation property holds true for Jh- 
(B2) For every w € Ji n H 2 , there exists an approximation r^w G Jh such that 

||w — r/iw|| + /i||V(w — rftW)|| < /i 2 1 1 W j 1 2 ■ 

The I? projection Ph : L 2 i-> satisfies the following properties (see [5]): for <f> G J^, 

(3.3) ||0 - P h <f>\\ + /»||VJ^0|| < Ch\\V<t>\\, 
and for G Ji n H 2 , 

(3.4) ||0 - P^H + ^||V(0 - P h( j>)\\ < Ch 2 \\A<f>\\. 

We now define the discrete operator A^ : *- >■ through the bilinear form a(-, •) as 

(3.5) a(v h , 4> h ) = (-A h v h , 4>) \/v h , 4> h G H h . 

Set the discrete analogue of the Stokes operator A = P{— A) as A^ = Ph{— Ah). Using Sobolev 
imbedding and Sobolev inequality, it is easy to prove the following Lemma 

Lemma 3.1. Suppose conditions (A.1 ), (131 ) and (B2 ) are satisfied. Then there exists a positive 
constant K such that for v, w, d) G Hh, the following holds: 



(3.6) 



(v • Vw,0)| < K < 



V2|| Vv ||V2||Vw||V2||A ft w||V2||0|| ) 
1 /2||A ft v|| 1 /2|| V w||||0||, 

1/2 ||Vv|| 1 /2||Vw||||0|| 1 /2||V0|| 1 /2 ) 

IIVwIIII^HA^II 1 / 2 , 

||Vw|| 1 /2|| A ^ w ||l/2||^ > ||l/2|| V ^ > ||l/2 



For examples of subspaces and Lh satisfying assumptions (Bl), (B2'), and (B2), we again 
would like to refer to [5] and references cited therein. 

Note that the semi-discrete solutions admit a priori and regularity estimates analogous to that 
of the continuous solution. Below, we present a Lemma containing certain estimates of the same 
and once again, we would like to refer to [5] for a proof. 

Lemma 3.2. Let < a < min{/xAi, 5}, where Ai > is the smallest eigenvalue of the Stokes' 
operator. Let the assumptions (A.\),(A!2),(R\) and (B2) hold. Then the semi-discrete Galerkin 
approximation U/j of the velocity u satisfies, for t > 0, 



(3.7) 

(3.8) 
(3.9) 



+ e 



-2at 



|Vu ft (t)|| + e 



-2at 



e 2as \\Vu h 



2 ds < K, 



e Zas \\A h u h (s)\\ z ds<K, 
(r^t^WAhUhmKK, 



where T*(t) = min{t, 1}. The positive constant K depends only on the given data. In particular, 
K is independent of h and t. 



The following semi-discrete error estimates are proved in [5]. 

Theorem 3.1. Let 0, be a convex polygon and let the conditions (A.1)-(A2) and (B1)-(B2) be 
satisfied. Further, let the discrete initial velocity u /i € Jh with u h = Ph^-o, where uq € Ji. 
Then, there exists a positive constant C , that depends only on the given data and the domain £1, 
such that for < T < oo with t E (0, T] 

||(u - u h )(t)\\ + h\\V(u - 0(t)|| < Ce ct hH- l l\ 



4 Nonlinear Galerkin Method 

In this section, we work with another space discretizing parameter H such that < h < H and 
both h, H tend to 0. We introduce two more spaces as follows: 

(4.1) 3 h = J H + J%, with 3 h = (I-P H )Jh 

Note that, by definition, the spaces Jh and J^ are orthogonal with respect to the L 2 -inner 
product (•, •). In practice, Jh corresponds to a coarse grid and J^ corresponds to a fine grid. 
The following properties are crucial for our error estimates. For a proof, we refer to pQ. 

(4-2) Hxll <cff|| X ||l, XG Jfc- 

And there exists < p < 1 independent of h and H such that 

(4.3) |o(0,x)l<(i-p)H0lli||Xhlli s 0eJ ff ,xeJf. 

From (I4.3|) . we can easily deduce that 

(4.4) p{H\\l + \\x\\i}<U + x h \\l ^Jff.xeJf. 

In the first modified nonlinear Galerkin method (NLG I), we look for a solution u h in Jh such 
that 

It is implemented on the interval (to,T], < to < T < oo. On the interval (0,to], classical 
Galerkin method is implemented, resulting in (uh,Ph)- 

Remark 4.1. Since the error analysis of NLG demands higher regularity of the unknown and 
this means higher singularity at t = 0, the idea is to avoid these kinds of singularity. 

For t > to, we look for (y H ,z h ) satisfying 

(yf , <f>) + pa(u h , cf>) + 6(u\ u h , 0) + / (3(t - s)a(u h (s), 0) ds = (f , 0), 

J to 

(4.5) M a(iA X ) + Ku h , u h , X )+ f Pit - s)a(u h (s), X ) ds = (f , X ), 

Jto 

for (f> e J H , X G Jf ■ We set y H (t ) = P H u h (to). 

In the second modified nonlinear Galerkin method (NLG II), we look for a solution u' 1 in 
such that 

u = y +z GJ^ + J/j 

satisfying 



(yf,0)+/ia(u fe ,^) + 6(u h ,u fe ,^)+ /" 



/3(t-s)a(u h ( S ),<£) ds = (f,0), 

(4.6) M a(ir\ X ) + 6(u h , x) + 6(y ff , z\ x) + /" /3(t - s)a(u h ( S ), x) ^ = (f , X ), 

Jto 

for G J H , x G . 



Remark 4.2. Here and henceforth, subscript means the classical Galerkin method and super- 
script means the nonlinear Galerkin method. 

Remark 4.3. Both the NLGs can be heuristically derived from as follows. We split the 

Galerkin approximation with the help of the L 2 projection Ph- 

(4.7) u h = P H u h + (I - P H )u h = y H + Zh, 

And we project the system hS. #)) on J#,J^ to obtain the coupled system: 

, 4g s f (yht,4>)+ ^a(u h ,(f)) + j^(3(t - s)a(u h (s),d>) ds = -b(u h ,u h ,4>) + (f,0), 

/or G J//, % ^ J?- Assuming the time derivative and higher space derivatives ofz h are small, 
various (modified) NLG methods are defined. In case, time derivative of z h is retained in the 
equation, different time-steps can be employed for the two equations, (much) smaller time-step 
for the equation involving y , for example kl = k, k2 = kl + 1, I £ N. In other words, z h 
remains steady as y H evolves for I times at each turn (see J2/), which in a way is equivalent to 
treating an evolution equation in y H and a steady equation in z h . 

The well-posedness of both the NLGs follows easily as in the case of Navier-Stokes equa- 
tions, see |15l [T]. But for the sake of completeness, we present below a priori estimates of the 
approximate solution pair {y^, z h }. And for the sake of brevity, we only sketch a proof (similar 
to the proofs in [5]). Note that the proof will go through for both the NLGs. 

Lemma 4.1. Under the assumptions of Lemma \3.2\ and for y H (to) = PH u h(to), the solution 
pair {y H ,z h } of (|^.5| ) (or \4-ff) ) satisfies, for t > to 

(4.9) ||y^|| 2 + e- 2Qi [ l e 2as \\Vu h {s)\\ 2 ds < K. 

Jt 

And if H is small enough to satisfy 

(4.10) fi-cL 2 H \\y H \\ 2 > 0, 
where Lh ~ [log h\ 1 / 2 , then the following estimate holds 

(4.11) \\z h \\<K. 

The constant K > depends only on the given data. In particular, K is independent of h,H 
and t. 

Proof. Choose <f) = Y H :X = z h in (|4.5p (or (I4.6P ) and add the resulting equations. Multiply 
by e 2at , integrate from to to t. Drop the resulting double integral term due to non- negativity. 
Multiply by e~ 2at to obtain fl4~9j) . 

For the second estimate, we choose x = z h in the second equation of (|4.5p (or ()4.6|) ) to 
obtain 

(4.12) ^||Vz ft || 2 <2^||Vy"|| 2 + — ^— f e- 2a ^\\Vu h (s)\\ 2 ds-b(y H +z\y H ,z h ). 

Use (|4.9p to bound the second term on the right hand-side of (I4.12p . Using (13. If) and (14. 21) . we 
find 

b(y H ,y H ,z h ) < 2 1 / 2 || y H || 1 / 2 || Vy H || 3 / 2 ||z / l 1 / 2 ||Vz /l || 1 / 2 < c\\y H \\\\Vy H \\\\Vz h \\ 
< ||| Vz h||2 + c || y H||2|| Vy ff||2 j 

b(z h ,y H ,z h ) < ||z ft ||||Vz h ||||y ff || 00 < cLH||z h ||||Vz h ||||Vy H || 
^^HVz^f + cL^llz^fllVy^ll 2 . 
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We have used the finite dimensional case of Brezis-Gallouet inequality (see |14l (3.12)]) 

lluftHoo < cL h ||Vu h ||; L h ~ [log h\ 1/2 . 
Therefore, we obtain from ()4.12|) 

|||Vz fe || 2 < K + 2^\\Vy H \\ 2 + c\\Vy H \\ 2 {\\y H \\ 2 + L 2 H \\z h \\ 2 ). 

Again with the help of f)3. 1 1) and ()4.2p . we note that 

H\\z h \\ 2 < KH 2 + K + e£ 2 j||y if || 2 ||z' l || 2 . 

Under the assumption on H, we have the required result, that is, (|4.1ip . This completes the 
proof. □ 

Remark 4.4. Further estimates of y H and z h can be obtained following the estimates proved 
in |5|/ for small enough H to satisfy 



5 Error Estimate 

In this section, we work out the error between classical Galerkin approximation and nonlinear 
Galerkin approximation of velocity. 

Before actually working out the error estimates, we present below the Lemma involving the 
estimate of z^. For a proof, we refer to 

Lemma 5.1. Under the assumptions Lemma 3.2 and for the solution of h3.2\) . the following 
estimates are satisfied for = (I — Pf/)u^ and for t > to 

( WzhW+HWzhW^K^H 2 , 
(5.1) \ \\z ht \\+H\\x ht \\ x <K(t)H 2 , 

( \\z htt \\ + H\\z htt \\i<K(t)H 2 . 

Remark 5.1. We note that Lemma \5. 1\ requires the estimates of higher order time derivatives 
of error due to Galerkin approximation: 

|| (u - u h ) (i) \\ + h\\(u - u h ) (i) ||i < K{t)h 2 , < i < 2, 

where (-)® means i th time derivative. For i = 0, the result is stated in Theorem \3.1[ For the 
remaining cases, the proof consists of differentiating the error equation in time and deriving 
estimates. Proofs are technical and lengthy and hence are skipped for the sake of brevity. 

In order to separate the effect of the nonlinearity in the error, we introduce 

u(€ J h ) = PffU + (/ - P H )u = y + z G 3 H + 

satisfying the following linearized system (t > to) 



(y t , 0)+/ia(u, <f>) + / P(t - s)a(u(s),<f>) ds = -b(u h , u h , + (f, <j>) <f> € 3 H 

(5.2) /ia(u, x) + / P(t - s)o(u(a), x) da = -b{u h , u h , X ) + (f , x) X € jf , 

J t 

and y(to) = yif(to)- Being linear it is easy to establish the well-posedness of the above system 
and the following estimates. 

Lemma 5.2. Under the assumptions of Lemma 3.2, we have 

\\Vu h \\ 2 + e~ 2at f e 2at \\A h u h \\ 2 ds < K 

Jto 



where the constant depends on uq and f . 



Ahu h \\ < K, 
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5.1 NLG I 

We define 

e := u h - u h = (y H - y H ) + (z h - z h ) =: e 1 + e 2 . 
We further split the errors as follows: 



(5.3) 



ei = y H - y H = {yh - y) - (y 11 - y) = £1 - Vi e J# 
e 2 = z h -z h = (z h - z) - (z h - z) = £ 2 - V2 e Jf ■ 



For the sake of simplicity, we write 

For the equations in £ and 77: subtract (15. 2p from (14. 8p and subtract (I5.2p from (I4.5P to obtain 

(5 4) f (£i, t , 0)+ 0) + / t * - a)a(£(a), 0) cfe = 

\ + / to /3(i- s)o(^(s),%) ds = -(z w ,x), 

^ 55 ^ f (Vi,u<l>)+ fJ-a(rj,(j))+J^(t-s)a(rj(s),^)ds = b(u h ,u h ,(l))-b(u h ,u h ,(j)) 
\ /-ia(v,x) + f to P{t-s)a{r](s),x) ds = b{u h ,u h ,x)-b(u h ,u h ,x), 

for G J# and x S J^. 

Lemma 5.3. Under the assumptions of Lemma 3.2 

(5.6) e- 2Qi / e 2aT \\i{r)\\ 2 dT <K(t)H 8 . 

Proof. Choose d> = e 2at £i, x = e 2crf £ 2 i n (15.4|) . add the two resulting equations and with the 
notation £ = e at £, we get 

(5 - 7) ijl 11 ^ 1112 ~ "H^ill 2 + ^II^N? + - s)e° i(£ - s) «(^(s),^) rfs < e^llz^H ||^ 2 || 

Using (|4.2p and (|5.ip . we can bound the right-hand side as: 

< e at .K(t)H 2 .cH\\£ 2 \\ 1 < ^-\\l 2 \\\ + K{t)H G .e 2at . 



And using (|4.4p . we have 



-a||^il| 2 + /^||^|| 2 >(w-f )\\ii\\l + w\\l2\\l 



As a result, we obtain from ([57 
d , : i , : o a 



i|| 2 + 2( W - f )||i|| 2 + wll^ll 2 + 2 /%(t - S )e^- s )a(^( S )i) ds < K{t)H & e 2at . 



'to 

Integrate from to to t and multiply the resulting inequality by e~ 2at . Note that the double 
integral turns out to be non-negative and hence 

(5.8) ||d|| 2 + e~ 2at fiU.Wl + \\L\\i) ds < K{t)H % . 

Jto 



Q 



To obtain L 2 (L 2 ) -norm estimate, we consider the following discrete backward problem: for fixed 

'h 



to, let w(t) G Jh, w = wi + W2 such that wi G Jh, W2 G be the unique solution of 



(t < t < t) 

( (0,w 1)T ) - /xa(^w)- f* P(s - T)a(<f>Ms)) ds = e 2a \<j>^ x ) 
(5-9) ^ - a(x,w)-/ r '/3( S -r)a( X ,w( S ))d S = e M ( X ^ 2 ) 

I wi(t) = 0. 

With change of variable, we can make it a forward problem and it turns out to be a linearized 
version of (|4,5p and hence is well-posed. As in [5], we can easily obtain the following regularity 
result. 

(5.10) I e- 2aT \\vr\\ 2 2 dT < C f \\k\\ 2 dr. 

J to J to 

Now, choose <fi = £ 1; x = £2 an d use (j5.4j) with <^ = wi, % = W2 to find that 
llt(r)|| 2 = (£i,w ljT ) -/xa(£,w) - |^( S -r)a(^w( S )) 

-~dt^ 1,W ^ + J t / 3 ( r - s ) a (^( s )' w ) ds -_/ /0(s-r)a(6w(fl)) ds + (z w ,w 2 ). 
Integrate from to to t and note that the double integral terms cancel each other. 

(5.H) / ||^(r)|| 2 dT = ((^(t),wi(t)) - «i(t ),wi(to)) + / (%, w 2 )dr. 

J to ^ *o 

But wi(t) = and £i(io) = yii(to) — y(*o) = 0. Next, we observe that 

wi G J // => Pf/Wi = wi, w 2 G => Ph w 2 = 0. 

As a result 

wi — PfjWy = 0, w 2 — i 3 f/w 2 = w 2 , or w — -Pj/w = w 2 . 

Therefore, 

(z ht , w 2 ) = (z ftt , w - PffW) 
(5.12) < ||z w ||||w -Pffw|| < if (t)H 2 .cH 2 1| w|| 2 . 

Prom (15. lip , we get 

f* - 1/2 



/ U(T)\\ 2 dT < K(t)e 2at H 4 ( [ e- 2ar \\w\\ 2 dT 

J to ^ J to 



I to J to 

Use (IBTTOl) to conclude. □ 

In order to obtain optimal L°°(L 2 ) estimate, we would like to introduce Stokes- Volterra type 
projections (Sh,S^) for t > to defined as below: 

Sh ■ Jh —> Jh, Sfr : Jh — > J^ , 

and with the notations 

Ci : = YH ~ S H u h G J H , C 2 := Zh - Sffuh G J^ 
the following system is satisfied. 

fia{C, 0) + /* P(t - s)a(C(s), 0) ds = 0, G J Hl 

For the sake of convenience, we have written C = Ci + C 2 - 



(5.13) 



MC, X) + It W ~ s)a(C(s), X ) ds = -(z ht , X ), X G Jf • 



1 n 



Lemma 5.4. Under the assumptions of Lemma 3.2 
(5-14) HCll + llCtll <K(t)H\ 

Proof. Choose <f> = e 2at Ci, X = e 2Q *C 2 to obtain 

(5-15) m||CI|?+ f /3(t-s)e a ^a(C(s)X)<e at \\z M \\\\C 2 \\. 

Jtn 



As in dSZZD-dSISD, we establish 

(5.16) e~ 2at f WCWjds < K(t)H 6 . 

Jtn 



'to 

Now from (|5.15|) . we have 

/*IICIIi < IICIIi I P{t - s)e a ^\\C(s)\\ 1 ds + |M|||C 2 ||. 

J to 

Use (|5.16|) to conclude that 

(5-17) HCIIi <K(t)H 3 . 

In order to obtain optimal Z°°(L 2 )-norm estimate, we would use Aubin-Nitsche duality argument. 
For that purpose, we consider the following Galerkin approximation of steady Stoke problem: 
let w/j G be the solution of 

(ia(v, w h ) = (v, Ci + C2)) v e J h- 
Writing wi = PffW/,, W2 = (I — Pf/)w/ l , we split the above equation as 



(5.18) 



fj,a(4>,-w h ) = (0,Ci), </> G Jff, 
Ma(%,w h ) = (x,( 2 ), XGJf. 



It is easy to establish the following regularity result. 

(5.19) \\wh = wi + W2II2 < c||Ci + £2!!- 

Now, put <fi = Ci, X = C2 i n (|5.18p and use (|5. 13[) with <f> = wi, x = w 2 to find that 
IICII 2 = Va(C,w h ) 

= -e at (z ht ,w 2 ) - f (3(t-s)e a ^a(t(s),w h ). 
Jt 

Use the fact that A^w^ = Ci + C 2 to obtain 

(5.20) IKII 2 + f P(t-s)e a(t - s \t(s)X) ds = -e at (x hu w 2 ). 

Jtn 



We can estimate the right-hand side as in f|5. 12|) . Integrate and observe that the double integral 
is non-negative to conclude that 



t 

(5.21) e~ 2at I \\tfds < K{t)H 



to 



1 1 



Now from (|5.20p . we find that 

||CH 2 <e Q *|(zM,w 2 )| + ||C|| I' f3(t-s)e a ^\\C\\ ds. 
Use (f5TT2D and (iOTT) to establish 



'to 



< K{t)H 4 . 

For the remaining part, we differentiate (|5. 13[) . 

f /ia(C t , 0) + /?(0)a(C, cf>) + j£ ft(t - s)a(C(s), <f>) ds = 0, G J H , 

1 ' J \ MC,x) + /3(0)a(C,x) + // /3 i (t-s)a(C(s),x) ds = -(z ha , X ), X € J^. 

Choose 4> = e 2at Ci,ti X = e 2Qi C 2 ,t to obtain 

^e 2a 1C t ||?</3(0)e at |K||i||Clli + ^ a 'IICIIi r^(t-«)e«( t -)||C(a)||icfe + e« t Ktt||||C !l , t ||. 

•'to 

Use kickback argument to find that 

llCtll? < cllCII? + ce~ 2at f \Us)\\\ds + K(t)H & - 
From (|5.16p - (|5.17p . we conclude 

(5.23) HCJi < K(t)H 3 . 

For L 2 -norm estimate, we again make use of Aubin-Nitsche duality argument. Similar to the 
problem (|5.18p . we consider 



(5.24) 



fia(4>,w h ) = (4>,e at Ci )t ), <f> e J#, 
fJ-a(x,^h) = (x,e at C2,t), X € Jf ■ 



It is easy to establish the following regularity result. 

(5.25) ||w ft = wi + w 2 || 2 < ce at \\C ht + C 2 ,tll- 

Now, put 4> = e Q *Ci, t , X = e at C2,t in (jSM]) and use (pQ2j) with <\> = e oi wi, x = e Ql w 2 to find 
that 

e 2at \\Ct\\ 2 = ne at a(C t ,w h ) 

= -e at (z htt ,w 2 ) - /3(0)a(C, w ft ) + «5 / /3(t - sK^a^s), w h ) 

<{e a iz, t J+c||C||+c||C||^e- (<5 - a)(t - s) d S }||w,|| 2 . 

Use (|5.25p and the estimate for ||£|| to establish 

K t \\<K(t)H\ 

□ 
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Now we are in a position to estimate L°°(L 2 )-norm of £, that is, of £± and £ 2 - Using the 
definitions of £j,Cj, i = 1,2, we write 

£1 = y# -y = (yff - £ffu h ) - (y-s#u ft ) =: Ci 

£ 2 = z/i - z = (zh - u h ) - (z - 5f u fe ) =: C 2 - 02- 
From (I5.4p and (I5.13p . we have 



(0i,t,0) + na(0,4>) + fl P(t-8)a(0(s), cf>)ds = (Ci >t ,0), ^Jff, 
M«(0, X) + Jt„ 0(< " s)a(0(s), X ) <fo = 0, * € Jf • 



Lemma 5.5. Under the assumptions of Lemma 3.2 

Hill < K{t)H A . 
Proof. Put = e 2at 0i, x = e 2o *0 2 in (1536]) to find 

(5-27) ^l|0il! 2 " «l|0il! 2 + M||0||! + J*P(t - s)e a ^a(0(s),9) ds < e at ||Ci, t ||||0i 

We recall that the spaces Jh and are orthogonal in L 2 -inner product. That is, 

for0G J H , xeJfc, O,x)=0. 

Hence 

||0i|| 2 < ||0i|| 2 + ||0 2 || 2 = ||0|| 2 < ||£|| 2 + IKII 2 , IICi.JI 2 < llCtll 2 - 

And 

-a||0i|| 2 + //||0||f = (fi - aX^W^Wj + fi\\e 2 \\l 
As a result, after integrating (|5.27p with respect to time from to to t, we obtain 



(5.28) ii^f + r(H0iii 2 +ii0 2 ii 2 ) d S < ( r e 2at ncii 2 ^) i/2 ( Aii^f+ncii 2 ) ^ v ■ 

J t y Jt ' v J t 



As usual we have dropped the resulting double integral as it is non-negative. We now use 
Lemmas 15.31 and 15.41 to conclude from (]5.28p that 



(5.29) H0J 2 + e~ 2at / e^fllfliHf + ||0 2 ||f) ds < K(t)H . 

J t 

We again choose x = e 2 "*0 2 in (|5.26p to find 

M ||0 2 || 2 = -//a(0i,0 2 ) - f p{t-s)e a ^a{e{s),e 2 ) ds 

Jto 



Using kickback argument, we obtain 

l|02||i<||0i||i + C (^(||0i|| 2 + ||0 2 || 2 ) ds^' 2 

Since 6\ G Jh, we use inverse inequality and (I5.29P to note that 

||0l||l < cH- l \\6\\ < K(t)H 3 . 

Hence, we conclude that 

||0 2 ||i <K(t)H\ 



i 2 



Now use (|4.2p to see that 

(5.30) ||0 2 || < K(t)H 4 . 

Combining f|5.29[) - (|5.30j) . we establish 

||0|| < K{t)H A . 

Use triangle inequality and the estimates of £ an d to complete the proof. □ 
We are now left with the estimate of r/, the error due to the nonlinearity. 

Lemma 5.6. Under the assumptions of Lemma 3.2 and that H is small enough to satisfy M-IO ) 
and 

lip - 2#||u|| 2 > 0, n - H(\\u\\ 2 + ||y H || 2 ) > 0, 

we have 

\\(u h -u h )(t)\\<K(t)H\ 
Proof. We choose 4> = e 2at r]i, x = e 2crf ?7 2 m fl5.5[) . 

(5-31) ^|ll^iH 2 + Mll^lli + ^(t ~ s)e a ^a(f,(s),f ) ) ds = e 2at A h ( Vl , V2 ), 

where 

and 

^h,i(Vi) = H u h, Ufc, r/i) - b(u h , u\ ?7 X ) 

= &(£ - v, Uh, vi) + K u h, £ - v, Vi) - K£ - v,t - v, vi) 
^,2(^2) = H u h, uh, V2) - H uh , u \ V2) 

= &(£ - 77, Ufc, %) + 6(u ft , ^ - r?, T7 2 ) - &(£ - T7, £ - r?, 

Therefore 

(5.32) &h(Vi, V2) = b (€ ~ V, u> »7) + K u h, £, »?)■ 

We estimate the nonlinear terms as follows: 

b(rj 2 ,u,rj) < ||»7 2 ||i||u|| 2 (||?7 1 || + ||»y 2 ||) < h2llill u ll2||??ill + -ff||u|| 2 ||T7 2 ||^. 
Therefore, for e, e% > 0, 

H(Vi,r] 2 ) <ehlli + eih2lli + c(e)(||u h ||2 + ||u||2)||^|| 2 

+ C (6,6 1 )||u||l||, 7l || 2 + J f/||u||2h2!l?. 

Now, from (I5.3ip . we find that 

^Wmf + 2//p(||T) 1 || 2 + ||r) 2 || 2 ) + tf t P{t ~ s)e a ^a(f ] (s),f t ) ds < 2e\\f ] \\l + 2e 1 ||r) 2 || 2 
(5-33) + c(e)(K|| 2 + ||u|| 2 )||£|| 2 + cC^eOHulllll^f + 2F||u|| 2 ||*) 2 || 2 . 
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We choose e = e\ = pp and assume that H small enough such that 

fip-2H\\u\\ 2 > 

to obtain after integration 

hill 2 + e~ 2at f\{\\m\\l + \\V2Wl)) ds < K(t)H s + K f \\ m {s)\\ 2 ds. 
J to •'to 

Apply Gronwall's lemma to establish L°°(L 2 )-norm estimate of rji. We note that 

||r/i||i < cH^Wr]^ < K(t)H 3 . 
For rj 2 , we again put x = e 2at r] 2 in (|5.5|) . 

(5.34) M ||ty? = +e 2at A h , 2 (r ]2 ) - pafo, f, 2 ) - f f3(t - s)e a ^a(f,(s), f, 2 ) ds. 

Recall that 

A h,2(f? 2 ) = Ki - v, u, V2) + K u h, £-vi,Vi) + b (€ - v, mi*!*)- 

And 

&(£ - vi,u,V2) + K u h,€ - viiVi) < (U\\ + ItoillXNh + HufcyiMli 

&(£ - ViiVuVi) < (Uh + hilWIMklMli 

K~V 2 ,u,r] 2 ) + b(-r] 2 ,r] 1 ,r] 2 ) < H(\\u\\ 2 + hilWH^IIi- 
Note that \\fli\\ < ||y || + ||y||- And under the assumption 

^-tf(||u|| 2 + ||y H || 2 ) >0 

we easily obtain that 

||»/ 2 ||i <K(t)H 3 

and hence 

\\r)2\\<cH\\V2\\i<K{t)H A . 
Now, triangle inequality completes the proof. □ 

5.2 NLG II 

In this subsection, we deal with the error estimate for NLG II. As earlier, we split the error in 
two, that is, e = — u h = £ — rf. The equations and hence the estimates regarding £ remain 
same and are optimal in nature. The equation in r\ reads as follows: 



(5.35) 



(li.t, 4>) + M??, <t>) + J* P(t - s)a(r](s), <p) ds = b(u h , u h , <p) - b(u h , u h , (p) 
Va{v, X) + J to - s)a(rj(s), x) ds = b(u h , u h ,x) ~ b(u h , u h , x) 

+b(z h ,z h , X ), 

for (p £ Jh and x S Jjf . 



Lemma 5.7. Under the assumptions of Lemma \5.6[ we have 

\\{M h -M h )[t)\\<K{t)H*. 
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Proof. We choose = e 2at ri l , x = e 2at r] 2 in f)5.35|) . 
(5.36) 

\j t t + + f t P(t - s)e^a(f,(s),f,) ds = e 2at {A h (rj 1 ,T] 2 ) + b(z h , z h , r, 2 )} , 

Since z h = Zh + r\ 2 — £ 2 , we have 

b(z h , z h ,rj 2 ) = b(z h + rj 2 - £ 2 ,z h - £ 2 , ?7 2 ). 



Now 



H€2i z h -£2>V2) < ll^2lll(ll z ft||l + ll&IIOIbtolll 

b(r) 2 ,z h - £ 2 ,r} 2 ) = b(r] 2 ,z,rj 2 ) < ci?||z|| 2 ||»7 2 [|f 
K z h,z h -$, 2 ,r] 2 ) < ||z^|| 1/2 ||z ft ||i /2 (||z^||i + ||^ 2 l!i)l! T 72l| 1/2 |l } 72lli /2 



1/2 1 1 it 1/2 |, n n , i|l/2|i II 1/2 . uj. n \ 

< cffllzfcllidlzfcHi + ||£ 2 ||i)ll»felli + cff 1 / 2 ||z h ||i||»y 2 ||i.cff 1 / 2 (||z h ||i + ||£ 2 ||i) 

< ^Wllll^lll + c - H 'll z /il|lli^2lllll r 72lll- 

As earlier, for small H, we obtain 

\\Vif + e~ 2at A(||^||? + ||r) 2 || 2 )) ds < K(t)H 8 + K{t)H 2 \\z h \\\ + K f \\ Vl (s)\\ 2 ds, 

Jtn Jtn 



which results in 



That is 



*?il| 2 + e- 2at A(M + \\v 2 \\ 2 )) ds < K(t)H e . 



\\ Vl \\<K(t)H 3 , 1177x11! <K{t)H 2 . 
As in the previous section, using only the second equation of (15.35j) we can easily conclude that 

\\rj 2 \\<K(t)H 3 , ||77 2 ||i < K(t)H 2 . 

□ 

Remark 5.2. The analysis reveals that the decrease in the order of convergence is due to the 
presence of b(zh, z^, x) i n the error equation. So, whereas in NLG I, we keep the nonlinearity 
in both the equations, in NLG II, the second equation is made linear in z h by dropping the term 
b(zh,Zh,x) an d which in turn appears in the error equation and is responsible for bringing down 
the rate of convergence in the above analysis. 

5.3 Improved Error Estimate 

In this section, we try to improve the rate of convergence, using the technique of [17]. The 
same technique is applicable for NSE also, but it is not straightforward, as the estimate of the 
function /(u) in their semi- linear problem does not always hold for our /(u) (which involves the 
non-linear term) and we have to be careful in order to obtain similar results. 
First, we note that the second equation of f)4.6j) can be written as 

(5.37) z h = $(y"), 



1fi 



where <3? : J# — > J^. Using this, we can write the equation in $(y#), for x £ J/^- 
Ha(y H + $(yn), X) + &(yn + $(yn ), yn , x) + &(yn > $(y# ), x) 

(5.38) +/ /8(t-s)a((y H + *(yH))(a),x)<fa = (fiX). 

■/to 

Lemma 5.8. Under the assumptions of Lemma 3.2 and that H is small enough to satisfy 

|- c tfK||i>o, 

we have 

(5.39) \\z h - $(y H )|| + ff||z ft - *(yfr)||i < K{t)H A . 

Proof. With the notation 3> e := — <&(y#) € J^, we have, by deducting (|5.38p from the second 
equation of fj4.8|) 



Aia($ e , X) + / - fl)o(* e (fl), X) <fa = - (zfc, t , x) - b(u h , u h , x) + b(y H , $(yn ), x) 
■/to 

(5.40) +b(y H + HyH),y H ,x)- 

Put x = ^"e to obtain 

Hl*e||l + / P{t~ s)a($ e (s), $ e ) ds = - (Zft,t, $e) - 6(*e, U fcl $ e ) 
•/to 

(5.41) -6($(y H ),$(yH),$e). 
Note that 

"(z M ,$e) < ||zMl|||^e|| < K(t)# 3 ||$e||l 
-6($e,U h ,$ e ) < C||u h ||i||$ e ||||$e||l < cH\\ U h \\ 1 || $ e ||f 

-6(*(yff),*(yH),* e ) = - $ e ,z h ,$ e ) 

< (iiz.ii^ii^np + ii* e ||^ii* e |i;/ a )i| Bfc || 1 ||*.ii^ii*.ii;/ a 

<A-(t)fr 3 ||*e||l+cff||z h ||i||* e |ll. 

Therefore, from (|5.4ip . we find 

(5.42) ^||*e|li+ [ P{t-s)a($ e (s),<f> e ) ds<K(t)H 6 + cfl-||u h ||i||* e ||i 
1 Jt 

We have used the fact that ||zfe||i < ||u^||i + ||y_ff||i < c||u.ft||i. And assuming H to be small 
enough to satisfy 

|-c#|K||i >0 

we establish after integrating (|5.42p 

/ ||* e ||ids < K(t)H 6 . 
Jto 

Use this result to estimate the integral term on (|5,42p to conclude 

||$e||l < K(t)H 3 . 

And hence 

||$e|| < Ci?||^ e ||l < K(t)H 4 . 

This completes the proof. □ 
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Lemma 5.9. Under the assumptions Lemma \5.8l we have 

(5.43) \\e 2 \\l<K{t)H & + c\\e 1 \\\ + cte- 2at f e 2tM ||ei(s)||?ds 

Jto 

(5.44) ||e 2 || 2 < K{t)H 8 + aff 2 ||ei|| 2 + cH 2 te~ 2at [ e 2as ||ei(s)||fcZs. 

Jto 



Proof. Recall that e 2 = z h - z h = (z h - <&(yn)) - {z h - 3>(y#)). With the notation $ e = 
z h — <3?(y#), we have e 2 = <l> e — The equation in <3? e can be obtained by deduction (|5.38|) 
from the second equation of (|4.6p . 



fia(^ - ei , x) + f M - s)a($ e (s) -e 1 , X )ds = - b(u\ y H , X )~ b(y H ,z h , X ) 

Jto 

(5.45) +b(y H ,Hy H ),X) + KYH + HyH):YH,X)- 
Put % = 3> e to obtain 

Mpli + f £(* ~ s)a(^ e {s), $ e ) ds = fia(e u $ e ) + /" (3{t - s)a{e 1 (s), $ e ) 

Jto Jto 

(5.46) +6(iA e l5 $ e ) + b(ei - $ e , y H , $ e ) + 6(ei, $(y H ), $ e ) 
Note that 

&(u ft , ei, $ e ) = 6(u h - ei - $ e + $ e , e 1; $ e ) 

< ||u ft || 2 ||ei||||* e ||i + ||ei|| 2 ||$ e ||i + ||*e||i||ei||i||$ e ||i + cff||ei ||i ||$ e || 2 
(for the last estimate, we have used H^Hl 4 < cH 1 ^ 2 ||$ e ||i) 

6(ei - $ e , y H ,$ e ) = 6(ej - $ e , u h - z h , $ e ) 

< ||ei||||u h || 2 ||$ e ||i + ||ei||i||z h ||i||$ e ||i + cF||y // ||i||^e|| 2 
b( ei ,<S>(y H ),<Z> e ) < ||ei||i(||$ e ||i + |K||i)||$ e ||i. 

Therefore, we find from (|5.46p 

(ji - cHWynMl^W 2 + 2 f p{t - s)a($ e (s),$ e ) ds < C || ei || 2 + c(l + tf 2 )|| ei || 2 

Jt 

(5.47) +ce~ 2at [ e 2as || ei (s)|| 2 ds. 
Assuming H small enough to satisfy 

H-cH\\y H \\i > 

we have, after integration 

(5.48) e- 2at e 2as \\<5> e \\ 2 < K(t)H e + ce~ 2at C e 2as ||ei(s)|| 2 cZs. 

Jt -'to 

Use (|5.48p to estimate the integral term in (|5.47p to obtain 

||$ e || 2 < K{t)H & + c\\e 1 \\\ + cte- 2at I e 2as \\ex{s)\\\ds. 

Jto 



And so 



||$ e || 2 < K(t)H 8 + cF 2 || ei || 2 + cH 2 te~ 2at [ e 2as || ei (s)|| 2 ds. 

Jto 

Using triangle inequality, we complete the proof. □ 
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Remark 5.3. Since the linearized error £ (that is, £1,^2) * s optimal in nature, so we obtain 
from [543\)-[544\> 



\\ V2 \\l<K{t)H & + c\\ Vl \\l + cte-^ f 

Jtn 



e** s \\ Vl (s)\\lds 

to 



\\r) 2 \\ 2 < K(t)H 8 + cH 2 ^ 2 + cH 2 te~ 2at [ e^'H^WHi***- 

Jto 



Use Lemma 5.7 to find that 

(5.49) ||i7 2 |li<*(*)# 6 + c||»/i||? 

(5.50) \\ V2 f<K(t)H 8 + cH 2 \\ Vl \\ 2 . 

Following [17] , we introduce the operator : J^, — > satisfying 

(5.51) a(v-i?fv, X ) = 0, V X G Jf. 
With the notations 

||v||fl = ||(7 - R%)v\\i, (v,w) fl = a((I - <)v, (/ - flf )w), 
we have, from Lemma 4.1 of |17| . 

( 5 - 52 ) c l|| v l|l < ll v l|fl < c 2||v||l, 

where ci,c 2 are positive constants independent of h,H. And similar to Lemmas 4.5 and 4.6 of 
[17] , we find for cj> £ Jh 



(yf,0)+My",0h = (f, (7 -7^)0)-/ P(t-s)a(u h ,(I-R*)<t>)ds 

(5.53) -6(u ft , u h , (/ - <)0) - &(z h , z h , i?f </>) 
(y^t,0) + M(yff,0)jJ = (f,(/-i2f)0)-/ p{t-s)a{M h ,{I-R%)<j>)ds 

Jto 

(5.54) -6(uh, Uft, (J — #f )0) + (u M , Tif 0). 
Now, for € J#, w e write the equation in ei = y# — y^ as 

(e lltj </>)+//(ei, <f>) R = - [ 0(t - s)a(ei, (/ - i?f )0) + (u M , i?f 0) 

+6(z h , z h , R%<j>) - 6(ei + e 2 , U/l , (7 - Tif )0) - b(u h , ej + e 2 , (7 - i#)0) 

(5.55) +6(e! +e 2 ,ei + e 2 , (7 - )0). 

Lemma 5.10. Under the assumptions Lemma \5.8l we have 

l|ei||i+ / Ktfds < K(i)i7 6 . 

7 in 
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Proof. Put = ei t m (|5.55p and observe that 

( u h,t,Rh e l,t) = -^_( u h,t,Rh e l) — ( u h,tt,Rh e l) 



= -(u h)t ,R% ei ) - ((I - PH)u tt ,R^ ei ) - ((uh - u) tt ,i^ei) 
< ^(ufc.t.ijfeO + A-^lleilli, 



- / /3(t- S )a( ei ( S ),(/-i?f)e lit ) ds = ^{- / P(t-s)a(e 1 (s),(I-R^)e 1 ) ds} 

J to Jtq 

+ /?(0)o(ej, (J - Rg)ei) - 5 f 0(t - s)a(ei(s), (I - R% )e a ) efa, 

-6(ei + e 2 , u ft ,(J - Rh)ex,t) ~ K u h,ei + e 2 , (I - Rh) e i,t) 

< c(||ei||i + ||e 2 ||i)||u/j2||ei, t ||, 

6(ei + e 2 , e!+e 2 , (J - i?f )e llt ) = b{e x + e 2 , u h - u ft , (/ - i?f )ei, t ) 

< c(||ei||i + ||e 2 ||i)(||u/j|| 2 + ||u h || 2 )||ei jt || 

b(z h ,z h , i?f ei, t ) = 6(z h - e 2 , z ft - e 2 , i?f ei, t ) = —b(z h ,z h , i?f ei) - 6(z M , z ft , ftf ei) 
- b(z h ,z hjt , i?f ei) + 6(z/j - e 2 , -e 2 , R^e 1)t ) + 6(e 2 , z h ,R^e ljt ) 
-b(zh,t,^h,Rh e l) ~ K z h,Zh,t,Rh e l) < cFllz^llillz^llilleilli 
b(e 2 ,z h ,R% ei , t ) < c||e 2 ||V2||e 2 ||i/ 2 (||z,|| 1 ||<e lit || 1 /2|| jR ^ eiit ||i/2 + || Zh ||i/2|| Zh ||V 2 || fl H eill 

< cff||e2||i||zft||i||ei )t ||i < c||e2||i||zft||i||ei )t ||. 

Here, we have used that — i??)ei ! t|| < ||eit|| + cJEf||ei t||i < c||ej And now we have 

l|ei,*|| 2 + ^l|ei||| < tf(t)ff 3 ||ei||i + +c||e 1 || 1 (||e 1 || 1 + f fi(t - a)||ei( a )||i) 

^{(u hj t, fif ei) - jf 0(i - s)a(e 1 (s), (I - R%)e x ) ds + b(z h ,z h , R% ei)} 
( 5 - 56 ) + c (ll e illi + ||e 2 ||i)||ei !t || +cfl r ||z h ||i||z/ l)t ||i||ei||i + c||e 2 ||i||z fe ||i||ei )t ||. 

Integrate (15.56p .use <|5.52[) and the fact that ei(io) = to find 



ei||f + I \\e lit \\ 2 ds<K(t)H 6 + c [ (|| ei || 2 + ||e 2 || 2 ) ds + (u htU R% ei) 
- / /3(t - s)a( ei (s), (/ - i?f )ei) + 6(z fc) z h) i?f e x ). 



'to 

As earlier, we estimate the last three terms to obtain 

(5.57) l|ei|| 2 + I \\eiA?ds < K(t)H 6 + c f (||ei|| 2 + ||e 2 || 2 ) ds. 

Jto Jto 

We note from (|5.49p and triangle inequality that 

||e 2 || 2 < K(t)H e + H^ll 2 < K(t)H 6 + || ei || 2 . 

Therefore 



l e i|li+ / We^fds < K(t)H 6 + c f ||ei|| 2 ds. 

Jto Jto 



I to Jto 

Use Gronwall's lemma to complete the rest of the proof 



9(1 



Remark 5.4. The Lemma \5.10\ tells us that 

\\r,x\\x<K(t)H 3 , 
and as a result, from Remark \5.3l we have 

\\r,2\\ + H\\ri2\\i<K(t)H 4 - 
Another application of triangle inequality results in 

||e 2 || + #||e 2 ||i < K(t)H A . 
For the final estimate, we write down the error equations in terms of e$, i = 1,2. 



(5.58) 



(ei, t , cj>) + na(e, </>) + L P(t - s)a(e(s), 4>) ds = -b(u h , u h , 4>) + b(u h , u h , (p) 



fia(e, x) + f to P(t - s)a{e(s),x) ds = -(z ht , x) - b(u h , u h , x) 



+ b(u h ,u h ,x)-b(z h ,z\x), 
Lemma 5.11. Under the assumptions Lemma \5.10\, we have 

\\(u h -u h )(t)\\<K(t)H\ t>t . 

Proof. With the notation A# = Ph{— A^), we choose (f> = Aj^eij in the first equation of (|5.58p 
to find 

u, d 

ll e i,t||-i + 2 " 6l " 2 + J ^ ~ • s )( e i( s )' e i.<) ds =b ( e i + V2,u h , A]j ei )t ) 
(5.59) + &(u / \e 1 + e 2 ,A^ 1 e 1 , i ). 

Observe that 

j t [ f* p(t - s)(ex(8),ex) ds} = f p(t - s)(ex{s), e ljt ) ds + /3(0)|| ei || 2 

-5 / j 9(i-a)(e 1 (a),ei) ds. 
J t 

Hence, we obtain from (|5.59|) 



/i d 



e i,dl-i+- — ll e ill + 



d 



2 '"' "~ 112 " J 7 /8(t-«)(ei(«),ei) ds} + 5 p(t - s)(ex(s), e x ) ds 

>t } J to 



If ' " 

2dt n " di 

(5.60) = P(0)\\ex || 2 + 6(ei + e 2 , u h , A^ei, t ) + 6(u ft , 6l + e 2 , A^e^). 

As earlier, we have 

6(ei + e 2 , u h , A^ei^) + 6(u ft , e x + e 2 , A^e^) 
<c(||ex|| + ||e 2 ||)(||u h || 2 + ||u h || 2 )||e lit ||_i. 

Integrate (I5.60|) and use the above estimate to find 

||ei|| 2 + f \\ex,t\\lx < c I (||ei|| 2 + ||e 2 || 2 ) ds < K(t)H 8 + c I \\ex\\ 2 ds. 

J t() J to J to 

Apply Gronwall's lemma to conclude 

l|ei|| 2 + f ||ei,t|| 2 x < K{t)H 8 . 
Combining with the Remark 5.3, we have the desired result. 



□ 
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Remark 5.5. It is clear from our above analysis is that the linearized error between NLG 
approximation and Galerkin approximation is of order H A in L? -norm. However, non-linearized 
part of the error may not always be of same order. For example, if the equation in z h contains 
only b(y H ,y H ,x), then the non-linearized part of the error (i.e. the equation inr}) will contain 
additional terms like b(y H ,z h ,x) an d b(z h , y H ,x) apart from the non-linear terms of the second 
equation of 115. 35]) . And with one of these terms, we believe, we can only manage H 5 order of 
convergence in I? -norm. 

Acknowledge: The author would like to thank CAPES for financial grant. 
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